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Recently, the deformed fi-Bose gas model based on so called //-deformed oscillators was 
0^ . proposed. For that model, the intercepts of r-particle momentum correlation functions (cor- 

relation functions at coinciding momenta of particles) were treated for r = 2, 3, within certain 
order of approximation in fi. In this work, we derive for /i-Bose gas model the exact expres- 
sions for r-particle correlation function intercepts, for all r, through Lerch transcendent and 
find their asymptotics (as functions of /t). For 2-, 3-particle intercepts and deformed distri- 
bution function (a£afc) the dependence on particles momentum is presented graphically. 
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^ 1 Introduction 

q \ The deformed /i-Bose gas model [1] is a particular implementation of such an efficient in nonlin- 
ear quantum physics tools as deformed oscillators (for their definition, diverse models and some 
applications, see e.g. [2]). The latter provide a generalization of the quantum harmonic oscillator 
and are often used in effective descriptions of various nontrivial properties of nonlinear systems 
and/or phenomena. Due to an effective representation in terms of deformed oscillator (s) we may 
! come to a more simple system whose behavior approximates the behavior of the initial one. Main 
properties of a deformed oscillator and its distinction from the standard quantum harmonic oscil- 
lator are characterized by the corresponding deformation structure function 4>(N) (see [2-4] for its 
notion). Depending on <f>(N) one distinguishes many types of deformed oscillators. Say, the one 
used in the elaboration of /x-Bose gas model is called //-oscillator and was introduced in [5]. Be- 
sides this rather exotic /i-deformed oscillator (for which, (f)(N) = ^ N ) there are such well-known 
g-deformed oscillators as Arik-Coon [6] and Biedenharn-Macfarlane [7, 8] ones, the p, g-deformed 
oscillator defined in [9] and also the g-deformed Tamm-Dancoff oscillator [10, 11]. Let us note 
that the latter four deformed oscillator models belong to the class of Fibonacci oscillators [12] (i.e. 
such that their structure function (f> n = <f>(ri) satisfies Fibonacci-like three-term recurrence relation 
with constant coefficients). Unlike, the //-oscillator represents the class of so-called quasi-Fibonacci 
deformed oscillators treated in [13] and that implies some peculiarities in concrete applications. 

Besides the /i-Bose gas model there exist many other deformed Bose gas models in which 
deformed oscillators were utilized. To list a few, the g-Bose gas model was studied in [14] (see also 
the references therein), p, g-Bose gas model was investigated in [15, 16] and in some other works. 
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Thermodynamical aspects of deformed Bose gases were studied e.g. in [17-19]. It is worth to 
mention that deformed oscillators were also effectively used for modelling a gas of (nonpointlike, 
composite) particles or quasiparticles [20-22], in the study of Bose-Einstein condensation [23] and 
of phonon spectrum in A He [24], for the description of interacting particles systems [25]. As it 
was shown in [26, 27], it is possible to realize composite bosons in terms of deformed oscillators 
(deformed bosons), and to relate the inter-component entanglement just with the parameter of 
deformation [28]. In the context of n-particle correlations, deformed oscillators were applied in [14- 
16], and the corresponding intercepts of correlation functions (along with their asymptotics) were 
obtained as functions of deformation parameter. 

In the preceding paper [1] the approximate expressions to certain order in /i were obtained for 
the intercepts of the 2nd and 3rd order momentum correlation functions in the deformed /i-Bose 
gas model, and for the corresponding asymptotics the exact expressions were found. It turns out 
that the results for the rth order intercepts can be obtained in exact form and this constitutes main 
subject of the present work. In fact, we derive the corresponding exact formulae given in terms 
of Lerch Phi-function (Lerch transcendent) for the general case of r-particle correlation intercept 
and also illustrate some of them graphically, with special attention to the particular cases r = 2 
and r = 3. One can hope that present results and further development /i-Bose gas model, say, 
concerning thermodynamical aspects will find effective applications to diverse physical systems. 

2 Derivation of the exact expressions for intercepts of the 
correlation functions in /i-Bose gas model 

Our treatment concerns deformed /i-Bose gas model as the secondary quantized system of specially 
deformed bosons. That means that each mode k in the deformed Bose gas is described by the 
deformed oscillator with the creation/annihilation operators a k , ak, and the number operator Ak. 
The defining relations for such a system of deformed oscillators (deformed bosons) are given as 



where the concept of structure function 4>(N) is used. As seen, we deal with the system of deformed 
oscillators which are mode independent i.e. the operators corresponding to deformed oscillators 
in different modes commute. 

Like in [1], and similarly to the corresponding nondeformed analogs [29], the concerned in- 
tercepts of the rth order momentum correlation function at a given momentum k are defined 
as 




(1) 
(2) 
(3) 



[Ak, ay] = 5 k ya 



k- 



[A k , ay] 



A w (k)= uwi/_ L 

(a k a k ) r 

Symbol (...) denotes usual statistical average for a system with Hamiltonian H: 




(4) 




(5) 
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where /3 = (ksT) 1 , ks is Boltzmann constant. Using (1) and (3) the formula (4) for the inter- 
cepts A^(k) can be rewritten in terms of the structure function <f>(N) as 



{) ~ PkF ' 



(6) 



where the notation [N] = [N]^ = <f)(N) for the structure function is used. Since we deal with the 
u-deformed Bose gas model [1], for (p(N) we take the structure function of u-deformed oscillator: 



M*0 = l N i 



N 



1 + fiN 



, n > 0. 



The intercepts (6) depend on the choice of Hamiltonian H. We choose the Hamiltonian in the 
form 

H = J2^N k (7) 

k 

and assume that deformed boson energy depends only on the absolute value k = |k| of the 
momentum. Then, using the Fock-like basis in which iVk|nk) = ft-kl^k), <K^V"k)|^k) = <K^k)|^k), 
the mode independence, and taking into account (5)- (7), for the r-particle correlation function 
intercept we obtain 



'1 _ p -/3Fui)\~r+l sr^oo n 



n—r+1 



H-/i(n-r+l) 



exp [-^n] 



1, cu=cuk=cu(k). (8) 



As seen, the intercept depends on u, on the momentum absolute value k and temperature T. To 
proceed, we expand the product in the nth summand of the numerator series (in front of exponent) 
into the sum of simple fractions: 



r-l 



III 



n — / 



+ ji{n — I) u r 



1 r ~ 1 1 11 r_1 



A (r V) 



+ u(n - Z) ' 



(9) 



The coefficients Af\fj) obey the recurrence relations: 



Z = 0,...,r-1; 



4 r+1) (") = -i-E 



^ Mr " 



(10) 



These relations allow to find the coefficients A^[\ to list a few: 

4%) = -i; 

4 2 V) = -i--, 4 2 V) = -i + -; 

/i 



4 3 V) = -i-f-^, 4V) = -i + 4 4V) 
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Then, in view of (9)-(10) the infinite sum in the numerator of (8) can be rewritten as 

OO ^ j oo 

■ ... • ; r exp [— BHum) = — exp [— BHum) + 

1 + un 1 + un-r + l ^ L J if ^ ^ v J 

n=0 n=0 

^ r—1 oo 

Likewise for the sum in the denominator of (8) we obtain 

oo ^ oo 1 °° 1 

exp \—Bftwn\ = — exp \— /3Hum) H — Aq (u) exp I — /?/kjnl . (12) 

^—J 1 + un L At ^—J L J a ^ 1 + un L J 

n=0 ^ ^ n=0 ^ n=0 ^ 

As seen, the second series in the r.h.s. of (11) and (12) are of similar form. Therefore in order to 
calculate this expression it is enough to perform summation of the generic series, that yields: 



_ _ 1 _ ~n 1 1 

E = E ^ - »•«/») - »*(«■*• 

n=0 n=0 

a = l-//Z, & = //, Z = l,...,r-1, z = e~ ptlu , 



where $ is the Lerch transcendent, see e.g. [30]. 

As result we arrive at the following expression for the thermal average {a\a,k) 

( < e l3huj _ J\-1 _ Q 

(a\a k ) = U(N k )) = { 1 i « ' „* , (13) 

Vfc \A t - 1 - / u" 2 (l-e- /3 ^)$(e- /3 ^,l, y u- 1 ), /i>0, V ; 

and the rth order averages (deformed analogs of r-particle momentum distribution) 

r-1 

((4) r (a k y) = 1 + I T 1 (1- e-^) A^VM^, 1, /i" 1 - 0, r > 2, (14) 

1=0 

which constitute one of our main results. The limit £t — > applied to (14) has a peculiarity, so we 
calculate {(al) r (cLk) r ) for £t = separately. Starting from its definition we arrive at 



((4) r K) r )U=o = (1 - e"^) ^ n(n - l)...(n - r + l)e" 



n=0 

Using Abel's identity for the latter sum we obtain: 



oo , 

(l-e - ^) Vn(n-l)...(n-r+l)e-^ n = (l-e - ^) lim i MiV-l)...(iV-r+l) 



e 



-P?iw(N+l) _ i 



7V->oo 

n=0 



I v J "' K ' e -/3Hui _ I 

1 —/3hu>(n+l) i ^ 00 

£rn(n-l)...(n-r+2) ^ ~ = (i_ e -/*-)^_ £ n(n-l)...(n-r+2)e* = 

/3r^ _ 1 (( a l) r " 1 K) r ~ 1 )U=o = • • • = p^; = jt; 
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Thus, 

((4)W)U=o= (e ,J^ 1)r - (15) 

From (13)-(14) and eq.(4) we obtain our next result - the intercept A^A;) of rth order correlation 
function, 

Ajft*) = (l + /i-^l - e-^^fW^^, 1, A*" 1 " «))• 

■(l + /*- 1 (l-e-^)4 1) 0*)*(e- /s ' ,w ,l ) /*- 1 ))" r -l ) r = 2,3,..., if /i > 0, (16) 

Using (15), in the no-deformation limit \x — > we recover: 

\ { ; ) (k)\, =0 = r\-l. (17) 

The obtained formula (16) presents the exact general expression for the intercepts under consid- 
eration as it covers all the orders r > 2. 

3 Series expansions in fi 

In the case of small deformation parameter /i it may be of interest to have an expansion of \ { ;\k) 
in /i. For that goal we expand the generic series present in (11) and (12) as follows: 



E tT7 _ n e- Q " = ^> a (0/A a = /3fu, Z = 0,...,r-1. (18 
The coefficients c s (Z) of the expansion are found as 

o s oo oo oo 

cM) = -r^-Y" ; K e~ an \ =e- al Y(-l) s (n-lYe- ain - l) = (-l) s e~ al V n s e" an = 

v y V 9a/ ^ V<W l-e- Q x'V dx/ x-1 



n=—l 



For the derivative in the last expression we perform the following transformation (m — I + 1): 



d \ s x m / d yx m - 1 f d y 1 / d y ^4 ^ / d \ s 1 
dxJ x — 1 V dx / x — 1 \ da;/ i-l V dx / ' V dx / x — 1 

j=0 

d \ s 1 



dx/ x — 1 

By induction it can be checked that the latter derivative may be presented as the finite sum: 



where the coefficients g J s satisfy the recurrence relation 

The solution of this recurrence relation is expressed through the Stirling numbers of the second 
kind {J}: 



■J+l 

In view of the latter, the coefficients c s (l) take the form 



C(Z) = S«T»* + e-"«"' ±fe^ + e^-l)- £ )w^- iy« 



j=l j=0 

Let us write out several first coefficients: 
co(0) = 1 + 

d(0) 



1 



(e^-1) (e^-1) 2 ' 
13 2 



cs(0) 



1 7 12 6 



r e( 9fioj — 1) (e l3tVjJ — l) 2 (e- 3 ^ — l) 3 (eP^ — l) 4 ' 



1 15 50 60 



c 5 (0) 



1 31 180 390 360 120 



( e /3hw_i^ ^ e f3hw_iy r e ptkj_i\3 ( e pnu_\Y ( e pnu_]\z> ( e pnw_\\t 



1 63 602 2100 3360 840 

c &\y) — T~MZ 7T + r „Rhui 1T2 ^~ f „Bha, TvT ^ TZMH, Tvt ^ TZMHi i\r ~^ 



Using the "initial" coefficients c s (0), all other coefficients can be determined by the formula 

i 

Cs (l) = e-W^j'ePtoi + e -^ Cs (0). 

3=1 

Then the sums (11) and (12), the latter being a particular case of the former, have the following 
Taylor expansions: 

oo 1 1 oo 

E = -(1 - e^)" 1 + ^fMc s (0)^. (20) 

The series (19) and (20) can be shown to be divergent for any \i > and, so, their practical use 
is doubtful. 
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4 Particular cases and comparison with results from other 
deformed Bose-gas models 

Consider two particular cases r = 2 and r = 3, earlier considered in [1] within a prescribed 
approximation. In these two cases, for the intercepts A^ and A^f (0 < \i < |) we obtain: 

Af = (1 - e-^)- 1 {(l - e-^)- 1 -(- + ^W^, 1, jT 1 ) - (- - -i)* (e"^, 



/i /X"/ \jJL jJL' 

I - e ~^Y l - /^(e"^, 1, /i^ 1 ))" 2 - 1, (2i; 



Af = (i- e H^-a{(i_ e H^)-i_(I + ^ + ^ 

- (--^ + ^3)$(e-^,l, / x- 1 -2)}- ((l-e-^)- 1 -^- 1 ^-^,!,^- 1 ))" 3 -!. (22) 

It is worth also to consider the following specially constructed [31] function r^\k) (also con- 
sidered in [1, 16]), useful in the experimental context: 

ri3)(k) _ A< 3 >(fc)-3A< 2 >(fc) 

It implies cancellation of unwanted distortions and provides improved purity. With account of (23) 
and the exact expressions (21)-(22) the explicit result for r^\k) readily follows merely by substi- 
tution (so we don't reproduce it here). 

In case when deformed bosons (here of /i-Bose gas) are used like in [1, 14, 16] to describe 
relativistic particles, the energy is specified as t\Wk = Vm 2 + k 2 . Then, the dependence of (a' k a^), 
\ff\k,T), \p(k,T) and r^\k,T) on the momentum k = |k| for the values /i = 0.1, 0.2 of 
deformation parameter, the temperatures being T = 120, 180 MeV, is such as shown in Figures 
1, 2, 3 and 4. 

From fig. 1 we see that the curves of deformed distributions lie below the usual (fi = 0) 
nondeformed Bose-Einstein one; with the increase of deformation parameter \x or decrease of tem- 
perature T the curves go lower and lower. For the momentum dependence of the intercepts A^ (k), 
\ < jf\k), and the function r^\k) (figures 2-4) we observe the similar behavior with respect to \i - 
larger deformation parameters \x correspond to lower curves, while lower temperatures correspond 
to curves lying higher. Besides, for all the plots in figures 2-4 the dependence on the momentum k 
shows asymptotical tending to the corresponding constant values for the intercepts (given by (27) 
below), with the large momentum asymptotics for (k) stemming from (23). 

The remaining part of this section is devoted to a comparison of the results obtained above 
with the analogous results presented in earlier papers on other deformed Bose gas models. The 
comparison concerns the correlation intercepts along with their asymptotics found within /i-Bose 
gas model [1], and those found in deformed Bose gas model of another type. For instance, in the 
case of p, g-Bose gas model the exact expression for the deformed analog of r-particle distribution 
and for the intercept of the rth order correlation function obtained in [15] respectively read (denote 
[r) P J = [r] P , q [r- l] p , g ...[l] p ,,): 

<(4r w> = Xrt^-JL (24) 
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Figure 1: Dependence of the deformed distribu- 
tion (<fi(N)) = (alak) on the momentum k, for de- 
formation parameter values /i = (pure Bose case) 
and /j, = 0.1, 0.2 (deformed case) and temperatures 
T = 120, 180 MeV 
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Figure 2: Dependence of the intercept A^(fc) on 
the momentum k, for deformation parameter values 
H = 0.1, 0.2 and temperatures T = 120, 180 MeV. 
Asymptotes are found from (27) at [i = 0.1, 0.2. 
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Figure 3: Dependence of the intercept A/f (fc) on 
the momentum fc, for deformation parameter values 
/i = 0.1, 0.2 and temperatures T = 120, 180 MeV. 
Asymptotes are found from (27) at /i = 0.1, 0.2. 
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Figure 4: Dependence of the correlation function 
r^ 3 (fc) from (23) on the momentum fc, for deforma- 
tion parameter values fj, = 0.1, 0.2 and temperatures 
T = 120, 180 MeV. Asymptotes stem from (23) af- 
ter the substitution of (27). 



((4)W) 



-pYie^ -q) r 



-1 



n 



j=0 



1. 



(25) 



With respect to the corresponding result (16) of the present paper, the expression (25), even though 
involving two parameters of deformation, has a simpler form (expressed in terms of elementary 
functions, not special ones), and this is presumably connected with the fact that p, g-oscillator 
belongs to the Fibonacci class whereas //-oscillator represents the class of quasi-Fibonacci [13] 
oscillators (therefore its treatment is more involved). 



S 



Now let us examine the asymptotic behavior. It is seen that the correlation intercepts tend to 
certain constant values at large momenta. These asymptotical values of X^\k,T) depend on /i 
only, and can be obtained directly from the definition (8) by taking the limit 

v-oo n _ _ n-r+1 -phwn r, ,-/3Kur , 

x(r) _ ,■ Z^n=0l+ M n •- l +At (n-r+l) e 1-1™ W G + - 1 _ 

A p, asympt ~ 11111 ' 1 ~ 11111 



= (l + //)>],! -1, [r] M ! = [rUr-l] M ...[l] M . (26) 

In the numerator, the dominating n = r term is retained whereas in the denominator the n = 1 
term dominates. For r = 2 and r = 3 this result is in complete agreement with the corresponding 
asymptotical values of the /i-Bose gas intercepts A^ and A^ earlier found in [1]: 

^asympt = { l + 1 1 ) [ 2 ] A J ~ 1 = ^ ~ ^ ^asympt = i 1 + [ 3 U ~ 1 = (1 + 2n){l + 3/i) ' ^ 

It is instructive to compare the asymptotical values (26) of the rth order correlation intercepts for 
/i-Bose gas model with the corresponding asymptotics in the case of p, g-Bose gas model [15, 16], 
that is 

^pq, asympt Mp,<T — 1" (28) 

The distinction, expressed on the /i-Bose gas side by the extra factor (l+/i) r , is connected with the 
fact that for p, g-oscillator we have [l] p , 9 = 1, whereas for //-oscillator [1]^ = 7^ 1. The exact 

expression for the asymptotics of the function rpf(k) can be obtained like in [1] by substitution 
of (27) into (23). 



5 Concluding remarks 

This work describes further steps with respect to [1] - here we extend those results and find them 
in precise form. More specifically, within /i-Bose gas model we have obtained the exact expres- 
sions (16), (17) for the intercepts of rth order correlation functions for any r > 2, as well as 
the explicit formulae (13)-(15) for (deformed and non-deformed) one- and r-particle distributions. 
Recall that in [1] the corresponding formulae were given only for the cases r = 2 and r = 3, and 
within certain approximation. Here the exact expressions for these particular cases (r = 2 and 
r = 3) are presented in (21) and (22) above. We have also obtained exact formula which gives the 
asymptotical values for the rth order correlation functions intercepts, as the corresponding func- 
tions of /t. Taylor expansions in the deformation parameter /t for one- and r-particle (deformed) 
distributions are obtained, see (19) and (20). Further analysis shows that they turn out to be 
divergent for any fi > 0. 

The dependence on the momentum k of thermal average (a k ak), of the second and third order 
intercepts Xp, (k) and Xp, (k), as well as of the function r^ (fc), see (23), was illustrated in the 
corresponding figures 1-4. The comparison of the obtained results with those of previous papers 

shows that the plots of Xp (k), Xp (k) and rp (k) are in qualitative agreement with earlier results. 

(2) (3) 

Moreover, the corresponding large momentum asymptotics of Xp (k) and Xp (k) completely agree 
with those found in [1]. If we compare the asymptotics (26) with the analogous asymptotics (28) 
found, also for arbitrary order r, within the p, g-deformation [15] we observe the following. In the 
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both cases the corresponding deformation of r-factorial does appear, but the result (26) differs 
from the formula (28), see also [15], by presence of the extra factor (1 + fi) r - this is caused by 
differing values of the deformation structure functions at n = 1. At last, few words concerning 
the further research on the relevant issues. For more comprehensive treatment of the /i-Bose gas 
model studied here and in [1] it is desirable to calculate some thermodynamic quantities (e.g. the 
critical temperature of condensation, the equation of state in the form of virial expansion) for 
/i-Bose gas, like this was done in [32] for the p, g-Bose gas model. 
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